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Abstract After some historical remarks we discuss different criteria of dynamical 
stability of stars, and properties of the critical states where dynamical 
stability is lost, leading to collapse with formation of the neutron star 
or a black hole. At the end some observational and theoretical problems 
related to quark stars are discussed. 
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1. Introduction 

It was first noted by S Chandrasekhar in 1931 that the mass of the 
compact star at very large densities, when electrons become relativis- 
tically degenerate cannot exceed a value, which is now called as Chan- 
drasekhar limit for white dwarfs (Chandrasekhar, 1931). After discov- 
ery of the neutron in 1932 the idea of the neutron star existence was 
developed. The matter in these stars is so dense, that neutrons become 
relativistically degenerate, again leading to the maximum mass of the 
neutron star. There are two important differences between white dwarfs 
and neutron stars, influencing the stability and mass limit. Neutron stars 
are about 6 orders of magnitude denser, than white dwarfs, making nec- 
essary to use general relativity (GR) in description of gravitation, instead 
of newtonian gravity. In GR the gravitation force is growing faster than 
1/r^, what reduce a stability and decrease the mass limit. On the other 
hand neutrons interaction become important at densities exceeding the 
nuclear one. The repulsion forces between nucleons at large densities 
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lead to additional stabilization against collapse, increasing the limiting 
mass of the neutron star. The influence of the second factor occurs to 
be more important, so at present the limiting mass of the neutron star 
is taken about two times larger than that of the white dwarf. 

Discovery of the quark structure of the matter led to suggestion of 
possible existence of quark stars, even more compact than neutron ones. 
In presence of indefiniteness concerning quark structure of the matter it 
is not possible now to make definite statements about possibility of exis- 
tence or nonexistence of stable quark stars, observational and theoretical 
investigations on this topic arc still in progress. 

In this review I first make a historical excursus into the problem, 
mentioning the results of the key works. Several criteria of stability 
are discussed, with the main attention to the static criteria, and ener- 
getic method, which permit to obtain the conclusion about the stability 
(sometimes approximately) in a most simple way. Critical states of com- 
pact stars at the boundary of the dynamic stability are considered, at 
which the star is becoming instable in the process of energy losses, and 
collapse begins leading to formation of the neutron star or a black hole. 
Physical processes leading to a loss of stability are discussed. At the end 
some observations and theoretical problems connected with quark stars 
are considered. 

2. Early development of the theory of compact 



From the theory of polytropic gravitating gas spheres with the equa- 
tion of state P = Kp^ it was known (Emden, 1907), that at 7 = 4/3 the 
equilibrium mass has a unique value 



where M3 = 2.01824 is a non-dimensional mass, corresponding to the 
polytropic index n = 3, 7 = 1-|-^. Such a star is in a neutral equilib- 
rium at any density. Chandrasekhar (1931) first noticed, that at large 
densities the electron degeneracy is becoming ultrarelativistic with the 
equation of state 



For some unknown reason Chandrasekhar (1931) took "/i equals the 
molecular weight, 2.5, for a fully ionized material", corresponding to 
K = 3.619 • 10^^ and had obtained the limiting mass from (1) equal 



stars: 1931-1965 




(1) 




(2) 



Dynamic stability of compact stars 



3 



to 1.822 • 10^3g= 0.91 Mq. Landau (1932) noticed, that for most stable 
nuclei: He^, C^^, O^^ etc. the value of /x, which is the number of nucleons 
on one electron, is equal to 2, and he obtained the realistic value of the 
limiting mass = IAMq. 

Soon after the neutron discovery (Chadwick, 24 Feb. 1932, letter to 
N.Bohr), "Landau improvised the concept of neutron star" in discussion 
with Bohr (Rosenfeld, 1974). The modern conception of the neutron 
star origin was due to Baade and Zwicky (1934), which appeared as 
a result of their studies of supernovae explosions. Their conclusion: 
"With all reserve we advance the view that supernovae represent the 
transitions from ordinary stars into neutron stars, which in their final 
stages consist of extremely closely packed neutrons", is an example of 
classic astrophysical predictions, finally proved only 35 years later after 
discovery of the pulsar in the center of the Crab nebula (Lovelace et 
al., 1968; Comella et al., 1969). Investigation of properties of matter at 
nuclear and larger densities expected in neutron stars had been started 
in papers of Hund (1936), Landau (1938) and Gamow (1938). 

The first neutron star model in GR was constructed by Oppenheimer 
and Volkoff (1939). They have derived equation of spherically symmetric 
stellar equilibrium in GR using metric in the form 

ds^ = -e" dt^ + dr^ + r^{de'^ + sin^ 9 d(p'^) (3) 
with equilibrium equations 

dP G{p + P/(P){m + ATrr^P/(^) 



dm 
dr 



dr r^{l — 2Gm/rc^) ' 

= 47rpr2, p = pQ(^l + ^y p{R) = 0, M = m{R). 



(4) 



An ideal degenerate neutron gas at zero temperature was considered 
with energy En, pressure P„, and rest mass density pQ as 

mic> , , 6.860 • 10^5 , , 

2A7T^hrp p 



247r2^^po Po 
Po = = 6.105 • lO'V 



where 
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/(y) = y(2y2_3)y^27I + 3smh-iy, (6) 

5(2/) =3y(2y2 + l)y^r,ri_3sinh-iy, y = 

They had obtained the following parameters of the neutron stars of dif- 
ferent masses: 

Table 1. Total mass of the neutron star model versus its radius, calculated by Op- 
penheimer and Volkoff (1939). 



Total mass in Solar units 


Radius in km 


0.30 


21.1 


0.60 


13.3 


0.71 


9.5 


0.64 


6.8 


0.34 


3.1 



Between white dwarfs and neutron star there is a neutronization and 
transition from normal matter with electrons and nuclei to superdense 
matter consisting of neutrons and other strongly interacting particles: 
mesons and adrons. The continuous curve M{pc) in the whole region 
from white dwarfs to neutron stars was first constructed by Wheeler 
(1958), where only neutrons have been considered in the superdense 
phase, like in Oppenheimer and Volkoff (1939). The results of Wheeler 
(1958) are given in Fig.l. The first neutron star model with a realistic 
equation of state, including mesons, hyperons and nuclear interaction, 
was constructed by Cameron (1959). The following equilibrium reactions 
between strongly interacting particles had been taken into account 

n + n ^ n + n + 27r° ^ A° + A° - 2(177) McV, 
p + j9^j9 + p + 27r°^S+ + S+- 2(251) MeV, 
A° + AO ^ A° + AO + vfO ^ S° + - 2(77) MeV, 

n + n-^ra + p + TT" -^^j + E" - 256 MeV, (7) 



n + A" ^ p + 7r" + A° ^ p + H~ - 204 MeV, 

AO + AO ^ AO + AO + 27r0 ^ ^o + ^o - 2(195) MeV. 

The nuclear interaction was considered according to Skryme (1959), see 
Fig.2. This equation of state is not perfect, because at large densities it 
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violates the casuality principle, according to which the sound speed can- 
not exceed the light speed (see Zeldovich, 1961). Neutron star models 
obtained by Cameron (1959) are represented in Fig. 3. The most im- 
portant result of this work was the indication that at realistic equation 
of state the limiting mass of the neutron star may exceed the Chan- 
drasekhar mass limit for white dwarfs. Therefore after a loss of stability 
a stellar core with mass exceeding the Chandrasekhar mass limit may 
stop its contraction forming a stable neutron star with enormous energy 
output, according to Baade and Zwicky (1934). In the model of Oppen- 
heimer and Volkoff (1939) the collapse would not stop, because of the 
low mass limit of the neutron star. 

Cameron (1959) had obtained that models at high densities had mo- 
notonously decreasing mass asymptotically tending to a constant value. 
That was probably the result of a crudeness of calculations, because the 
behavior of the curve M{pc) at large densities has oscillating, and not 
monotonous dependence. Such oscillations had been present in calcula- 
tions of Ambartsumyan and Saakyan (1961), who also improved physical 
description of superdense matter (Ambartsumyan and Saakyan, 1960). 
The nature of high density oscillations of mass in GR equilibrium was 
explained by Dmitriyev and Kholin (1963) (sec also Misner and Zapol- 
sky, 1964). They had shown that all models beyond the first neutron 
star mass maximum on the curve M(pc) are unstable, and after each 
new extremum the model acquires one new unstable mode. Later this 
problem was analyzed at length in the book of Harrison et al. (1965). 
Detailed calculations of neutron star models with a realistic equation of 
state had been done by Saakyan and Vartanyan (1964). 

3. Criteria of hydrodynamic stability 

The exact approach to the problem of dynamic (linear) stability is 
based on the solution of the equations for small perturbations, and find- 
ing eigenvalues and eigenfunctions of these equations. In the conserva- 
tive system the variation principle may be derived, which determines 
the exact value of eigcnfrequency if the exact eigenfunction is known. 
In practice even use of an approximate linear eigenfunction often define 
the eigcnfrequency with a good precision. The variation principle for 
a spherically symmetric star in GR was first derived by Chandrasekhar 
(1964) from the equations of small perturbations, and later by Harrison 
et al. (1965) by variation of the potential energy of the star. The exact 
method for determination of the dynamic stability of the star to the 
adiabatic perturbations based on the sequence of static solutions was 
formulated by Zeldovich (1963). He had shown, that in the extremum of 
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MASS IN SOLAR UNITS 



Figure 1. The relation between central density and the mass of various degenerate 
star models. Chandrasekhar's curve is for white dwarfs with a mean molecular weight 

2 of atomic mass units. Rudkjobing's curve is the same except for inclusion of the 
relativistic spin-orbit effects Rudkjobing (1952). The curve labeled "Oppenheimer 
and VolkofP' is for a set of neutron star models. The solid line marked "Wheeler" is 
a set of models computed with a generalized equation of state, from Cameron (1959) 



the function M(pc), or M{R) along the sequence with a constant entropy 
S adds or subtract one unstable mode. On the curve of cold stars on 
Fig.l the stability is lost in the first maximum due to neutronization, or 
GR effects, but it is restored in the minimum, where stable neutron star 
appears. The static criterium of stability was generalized for the case 
of arbitrary rotating stars with arbitrary distribution of entropy over 
the star in the paper of Bisnovatyi-Kogan and Blinnikov (1974). It was 
noted, that in presence of three quantities conserved during adiabatic os- 
cillations, the critical points are the extrema of three functions M(pc)\s,j, 
s{pc)\M,j, J{Pc)\M,s. Here along the sequences not only total values of 
stellar angular momentum J and entropy S should be conserved, but 
also their distributions over the lagrangian coordinate N, which is the 
number of baryons inside the cylindrical radius, j{N) and s{N). Note 
that in the newtonian case lagrangian coordinate is coincides with la- 
grangian mass m = pN = poN, but differs from it in GR case. For 
polytropic equation of state and angular momentum distribution similar 
to rigid rotation 
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Figure 2. Skryme's equation of state. The hard core modification is shown at the 
upper right, and the dashed hne is the equation of state for a non-relativistic Fermi 
gas of neutrons, from Cameron (1959). 



P = Kp'+'/-, K = const • e^ j{m) = ^^[1 - (1 - mf'] (8) 

examples of application of the static criteria are represented in Fig. 4 
from Bisnovatyi-Kogan and Blinnikov (1974). 

Let us stress, that along the curve appropriate for a static criterium of 
stability the rotation does not remain rigid, but only angular momentum 
of each lagrangian loop should be conserved. It means that if you have a 
curve of rigidly rotating stars with a fixed total angular momentum you 
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MASS IM SOLAR UNITS 

Figure 3. The observable and proper masses of neutron star models with non-ideal 
matter, from Cameron (1959). 



should construct an adjusted curve from each point of this curve with 
a constant angular momentum of a lagrangian loop, which maximum 
indicated a loss of stability. The rigidly rotating model is on the border 
of stability, when coincides with the maximum of the adjusted curve. 
For a cold white dwarf with the equation of state 



p = Ax^ (2 + aix + a2x'^ + a^x^^ 



P = B 



X {2x'^ - 3) (x^ + 1) + 3 In (x + VTTx2) 



A = 9.82 X 10° g cm""*, B = 6.01 x 10^^ erg cm"''. 



,-3 



(9) 



-5 



1.376 X 10"^ 



ai = 1.255 X 10"^ a2 = 1.755 x 10 03 

the adiabatic index is close to n = 3 at the critical point, and the 
eigenf unction is close the the linear one. At homologous transformation 



Dynamic stability of compact stars 



9 



M 



K 




1.000 - 




2.10 - 



0.99S 



2.05 L 



log Pc 



J I L 

3 -2 -1 



-3 



-2 



-1 







-3 







Figure 4- The dependence of the mass M on the central density for fixed entropy and 
angular momentum, for isentropic models having a polytropic index n = 3.03. The 
angular momentum distribution is given in (8) (left). The dependence of K on central 
density for fixed angular momentum distribution, corresponding to the mininmm of 
mass in the left figure, for isentropic models with polytropic index n = 3.03 (right). 
Prom Bisnovatyi-Kogan and Blinnikov (1974). 

the rigid rotation is preserved, so the adjusted curve for white dwarfs 
differs only sHghtly from the curve of rigidly rotating models, see Fig. 5 
(left). For a more complicated equation of state with highly variable 
adiabatic index, like in the equation of state 



the value 71 = 4^ at p = pi and falls off steeply when p > pi. The results 
of calculations for a model with solid-body rotation and three associated 
series of models are given in Fig.5 (right). The loss of stability occurs at 
the point of intersection of the associated curve D with the solid-body 
one, this point coinciding with the maximum on the curve D. It can be 
seen from Fig.5 (right) that the point of stability loss differs from the 
maximum on the solid-body curve by almost 5% with respect to pc- An 
extension of static criteria to general relativity and toroidal magnetic 
fields is made in the paper of Bisnovatyi-Kogan and Blinnikov (1974). 
Static criteria remains valid also in the presence of phase transitions. 
Variational principle and other methods of investigation of stability in 
presence of phase transitions are considered by Bisnovatyi-Kogan, Blin- 
nikov and Shnol (1975). The review of stellar stability methods is given 
in the book of Bisnovatyi-Kogan (2002) . 
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Figure 5. The dependence of the mass M on the central density for rigidly (Curve 
A) and differentially (Curve B) rotating white dwarfs with the equation of state from 
(9). The angular momentum distribution of the Curve B is that of the extremal 
model of Curve A (left). The dependence of the mass M on the central density for 
rigidly (Curve A) and differentially (Curves B,C and D) rotating models with the 
equation of state from (10). The angular momentum distribution of the Curves B,C 
and D are those of the rigidly rotating models of the Curve A at the points of its 
intersection with the Curves B,C and D. The curve C is intersecting the Curve A in 
its maximum. The maximum of the Curve D coincides with the point of intersection, 
which is a critical point for rigidly rotating models. (right). Prom Bisnovatyi-Kogan 
and Blinnikov (1974). 



4. Energetic method 

In the case of a complicated equation of state and entropy distribution 
along the star, the static criteria is hard to apply. In this case it is 
more convenient to use an approximate energetic method, which follows 
from the exact variation principle where linear trial function is used for 
estimation of the eigenfrequency. 

Energetic method was first applied for investigation of stability of 
supermassive isentropic stars by Zeldovich and Novikov (1965) where 
first GR corrections to the energy of such star have been found. This 
method was generalized for arbitrary stars by Bisnovatyi-Kogan (1966). 
Equations have been derived, giving approximate description of equilib- 
rium and determining the critical point where stability is lost. The first 
equation represents equilibrium condition, following from first variation 
of the total potential energy of that star (gravitational plus internal for 
nonrotating stars) equal to zero. The model with adiabatic index n = 3, 
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7 = 4/3 is in a neutral stability with a linear eigenfunction. In appli- 
cation of the energetic method the density distribution of the model is 
prescribed by the polytropic n = 3 distribution, p = pc<f{m/M), where 
m is lagrangian mass coordinate, and ip{m/M) is connected with the 
Emden function corresponding to n = 3, with a linear trial function 
5r = ar, Q;=const. The equation of state is arbitrary with the pre- 
scribed entropy distribution over the lagrangian coordinate, and first 
GR correction to the energy, including all terms ~ R/Rg, R is stellar 
radius, Rg = 2GM/(? is the gravitational radius, are taken into account. 
The equilibrium equation in the approximate energetic method is: 



3p4/3 f pj^ _ 0.639 GMV3 _ 1 84 ^ ^ ^^^^ 

Jo V (m) ^ 

Under the above mentioned conditions the second variation of the en- 
ergy is also reduced to the integral relation. Zero value of this relation 
approximately corresponds to the critical state of the loss of stability. 
We have 

4^ *^L__,.84^l!^=0. (12) 

'^^ Jo \ 37 ip{m/M) c2 ^ ^ 

These equations had been obtained by Bisnovatyi-Kogan (1966), using 
the expression of the energy with the prescribed distributions of the 
density (Emden polytrope n = 3) and entropy (arbitrary) over the la- 
grangian coordinate m/M, 

' = /„ •^(O'^irn-l — 5.06-^ (13) 

= / £;(p, T) dm- 0.639 GM^/^ _ 0.918 

Jo c 

dm = Airpr'^dr 

by differentiation over the central density. Critical parameters of iron 
isentropic stellar cores for masses between 5 and 1000 Mq at the point 
of the loss of stability had been calculated by Bisnovatyi-Kogan and 
Kazhdan (1966). The dependence of the critical central density on mass 
is represented in Fig.6 (Bisnovatyi-Kogan, 2002) for a wide region of 
masses. At low masses, high densities the stability in the iron core is 
lost due to neutronization. For larger masses the temperature in the 
critical point is increasing, and stability is lost due to iron disintegra- 
tion, decreasing the adiabatic power, and making it 7 < 1. At masses 
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around 1000 Mq with decreasing density and temperature, the adiabatic 
power becomes less than 4/3 mainly due to pair creation, and at farther 
increasing of mass the effects of GR are the main reason of the loss of 
stability (Zeldovich and Novikov, 1965). At large masses the dependence 
Pc{M) is described by the relation 

Pe,cr(M)= 2.4 X 10^^-3 (^^j g cm'^ (14) 

where /j, is the molecular weight. 

Models of hot isentropic neutron stars had been calculated by Bisnova- 
tyi-Kogan (1968), where equilibrium between iron, protons and neutrons 
was calculated, and the ratio of protons and neutrons was taken in the 
approximation of zero chemical potential of neutrino. The stability was 
checked using a variational principle in full GR (Chandrasekhar, 1964) 
with a linear trial function. The results of calculations, showing the 
stability region of hot "neutron" stars are given in Fig. 7. Such stars 
may be called "neutron" only by convention, because they consist mainly 
from nucleons with almost equal presence of neutrons and protons. The 
maximum of the mass is about TOMq, but from comparison of the total 
energies of hot neutron stars with presupernova cores we may conclude, 
that only collapsing cores with masses less that 15Mq have a chance 
to stop at the state of a hot neutron star, and larger masses collapse 
directly to black holes. 

The energetic method was generalized for rotating stars in the paper of 
Bisnovatyi-Kogan, Zeldovich and Novikov (1967), and in a more precise 
formulation by Bisnovatyi-Kogan and Ruzmaikin (1973). The density 
distribution remains the same, and the rotational energy term, written 
as a function of the total angular momentum was added. Rotational 
energy of the whole star depends on the central density as ~ P^^^, 
like the first GR correction term in (13). Therefore next order GR 
corrections had been taken into account ~ R~^, which include also GR 
corrections to the internal and rotational energies. Calculations of the 
second GR corrections to the energy of nonrotating stars have been first 
done by Vartanyan (1972). Energetic method with two GR corrections 
for rotating stars was applied for investigation of stability of supermasive 
stars by Bisnovatyi-Kogan and Ruzmaikin (1973). It was obtained, that 
for large angular momenta rotation stabilizes the star against a collapse 
and prevent formation of a black hole. It is in correspondence with a 
fact, that Kerr solution for a rotating black hole exist only for sufficiently 
small angular momenta, and there are no black holes at larger values. 
Only radiation with 7 = 4/3 was taken into account in the equation of 
state, what lead to absence of equilibrium configurations at low angular 
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Figure 6. Central density as a function of stellar core mass for the critical states, 
from Bisnovatyi-Kogan (2002). 



momenta j. The curves for equilibrium configurations at different j 
are represented in Fig. 8 on the plane: equilibrium entropy Seq - x = 
Rg/R. The results are expressed in non-dimensional variables and are 
valid for an arbitrary mass star with the adiabatic equation of state 
P = Kp^l-K The energetic method gives asymptotically exact results 
for slowly rotating supermassive stars, but even at large j the results 
remain correct qualitatively.^ 

5. Search of quark stars 

At very high densities the equilibrium lowest energy state corresponds 
to quarks. However because of big indefiniteness of the theory the con- 
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Figure 7. Equilibrium stable configurations on the mass M, central rest-mass den- 
sity pc diagram for superdense isentropic stars (hatched). The stars become stable at 
sufficiently large densities because of nucleon pressure, and become unstable again at 
larger densities due to GR effects, from Bisnovatyi-Kogan (1968). 

dition at which the quark state becomes preferable lays in a big density 
interval. Moreover, the process of transition between nucleon and quark 
state is also not definite, and even at very large densities the time of 
this transition may be enormously long. First the idea of a possibility 
of quark stars was suggested by Ivanenko and Kurdgelaidze (1965). De- 
pending on the theory parameters there is a possibility of 3 types of quark 
stars: strange stars (SS), which consists of strange matter which may be 
stable at nuclear density with a zero pressure; hybrid stars with hadronic 
envelope, which has a core of pure quarks (QC), or has a mixed core of 
hadrons and quarks (MC). There are different approaches to describe a 
quark matter, see Aguirre and De Paoli (2002), Andersen and Strickland 
(2002), Burgio et al. (2002), Sedrakian and Blaschke (2002), Kohri et 
al. (2003), Lugones and Horvath (2003), Alford (2003), Mishustin et 
al. (2003). The results of calculation of hadronic (H) and quark stellar 
models (SS, QC and MC) in Hard-Dense-Loop approach are represented 
in Fig. 9 from Thoma et al. (2003), where one of a model parameters is 
changing. The free quarks exist in the state of deconfined quarks, and 
the density when deconfined quarks become energetically preferable is 
also rather indefinite (Berezhiani et al., 2003). 
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Figure 8. The curves Se.q which characterizes the equilibrium entropy versus of 
X = for different values of the angular momentum j. The broken part corresponds 
to the nonphysical states Seq < 0. The parts of curves to the left of minima correspond 
to the stable states, from Bisnovatyi-Kogan and Ruzmaikin (1973). 

Larger scattering of the properties of quark stars is obtained by An- 
dersen and Strickland (2002), who used the same Hard-Dense-Loop ap- 
proach, but with wider variation of parameters. Their models of quark 
stars are represented in Fig. 10. 

In presence of such problems with theory, only observations could give 
an answer about the possibility os existence of quark stars. Even in pres- 
ence of theoretical indefiniteness, it is impossible to have a neutron (pure 
hadronic) star with a radius smaller than about 10 kilometers. Observa- 
tional discovery of a compact star with considerably smaller radius would 
be an important evidence that this star is not pure hadronic. Such claim 
was made by Drake et al. (2002) basing on deep Chandra observations 
of the isolated compact star RX J1856. 5-3754. They wrote: "We argue 
that the derived interstellar medium neutral hydrogen column density 
of 8 X 10^^ < Nh < 1.1 X 10^0 cm favours the larger distance from 
two recent HST parallax analyses, placing RX J1856. 5-3754 at ~ 140 pc 
instead of ^ 60 pc. ... The combined observational evidence — a lack 



16 




8 10 12 14 

R[km] 



Figure 9. Mass-radius relation for pure strange quark matter stars (left) and hybrid 
stars (right). GO - G4 models of hybrid stars corresponding to different parameters 
of the model. H: pure hadron star, QC: star has a quark core, MC: star has a mixed 
core, from Thoma et al. (2003). 

of spectral and temporal features, and an implied radius i?oo = 3.8- 
8.2 km that is too small for current neutron star models — points to a 
more compact object, such as allowed for quark matter equations of 
state." Subsequent analysis of observational data did not confirm this 
conclusion. Thoma et al. (2003) came to contrary result of an unusu- 
ally large value of the radius even for a neutron star: "Combining such 
a X-ray spectrum with the optical spectrum, one finds a black-body 
radius Rod > 17 km, indicating a stiff equation of state, which would 
exclude a strange quark matter star or even hybrid star. This is a rather 
conservative lower limit for the radius since a black-body emitter is the 
most efficient radiator." Similar conclusion was done by Turolla et al. 
(2003): "While a quark star may be a conceivable option, present obser- 
vations do not necessarily demand this solution and more conventional 
scenarios involving a neutron star are certainly possible. Neutron star 
models based on a two-temperature surface distribution can account for 
both the X-ray and optical emission, giving at the same time acceptable 
values for the stellar radius." Meanwhile a new candidate for a quark 
(strange) star in the X-ray source 4U 1728-34 was suggested by Bombaci 
(2003), and farther discussion about this object is expected. 

So far no reliable candidate for a quark (strange) star was found by 
observations. The theoretical predictions remain uncertain, and could 
explain either existence or nonexistence of quark stars. The solution 
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Figure 10. Mass-radius relation for a quark star with A//i = 1.6 and A//i = 1. 
The weak-coupUng results for the same choice of renormalization scales are shown as 
dashed lines, from Andersen and Strickland (2002). 

of intriguing problem of cosmological gamma ray burst (GRB) origin 
was proposed by Berezhiani et al (2003) basing on energy production 
during transformation of the neutron star into a stable quark (strange) 
star. This model could explain the connection of GRB with supernovae 
explosions with a formation of a neutron star, and subsequent huge 
energy release during transition to a state of a quark star producing 
GRB. An attractive feature of this model is a possibility of explanation 
of any time delay between SN and GRB explosions, which depends on 
the transition time from hadron to quark state, and is very sensitive 
to parameters. In this model the time delay may be arbitrary long, 
so majority of SN explosions in accordance with observations do not 
produce GRB, because this time exceeds the cosmological Hubble time. 

6. Summary 

1. Formation of quark (strange) stars does not follow unambigu- 
ously from the theory, which may be compatible either with existence 
or nonexistence of these objects. 
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2. There are no firm observational contradictions to the conventional 
model of the hadronic neutron star. 

3. Speculations about connection of the cosmological GRB with tran- 
sition from the hadronic to quark star seems to be interesting, because 
they explain connection between GRB and supernovae explosion, with 
arbitrary time delay between these events, including very large, exceed- 
ing the Hubble time. 
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Notes 

1. Because an approximate trial function does not give an absolute minimum of the energy 
functional, the loss of stability happens, strictly speaking, "before" the point obtained by the 
energetic method (EM), when going along the sequence of decreasing entropy or angular 
momentum. It means, that the model is definitely unstable in the EM critical point. 
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